Field dependence of the vortex structure in d-wave and s-wave superconductors 
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We study the vortex structure and its field dependence within the framework of the quasi-classical 
Eilenberger theory to find the difference between the d^2_y2- and s-wave pairings. We clarify the 
effect of the d^2 _y2-wave nature and the vortex lattice effect on the vortex structure of the pair 
potential, the internal field and the local density of states. The d^2_y2-wave pairing introduces 
a fourfold-symmetric structure around each vortex core. With increasing field, their contribution 
becomes significant to the whole structure of the vortex lattice state, depending on the vortex 
lattice's configuration. It is reflected in the form factor of the internal field, which may be detected 
by small angle neutron scattering, or the resonance line shape of fiSK and NMR experiments. We also 
study the induced s- and d^^y-wave components around the vortex in d^2_y2-'wa,Ye superconductors. 

PACS numbers: 74.60.Ec, 74.60. Ge, 74.72.-h, 76.75.+i. 



I. INTRODUCTION 

Much attention has been paid to the vortex structure 
in high- Tc superconductors. Many researchers try to de- 
tect the dx2-y2-wave nature of the superconductivity in 
the vortex structure of the mixed state. It is necessary 
to clarify how the difference between the d^2_y2-wave 
and conventional s-wave pairings appears in the mixed 
state. The point is how the vortex structure is affected 
by the anisotropy of the energy gap in the fia;2_j,2-wave 
pairing, particularly by its nodal structure. It is known 
that the difference appears in the field (H) dependence of 
the zero-energy density of states (DOS) N{0). Voloviktl 
theoretically suggested that the dependence is given by 
A^(0) oc Vi? in the da:2_y2-weive pairing and A^(0) oc H 
in the s-wave pairing. The difference between the two 
pairings was indeed confirmed by the self-consistent cal- 
culation of the quasi-classical Eilenberger theory for the 
vortex lattice case, but the field dependence of ^(0) de- 
viates from the exact or i?-linear behavior .□ In this 
paper, by using the quasi-classical method, we discuss 
other dx2-y2-'wave natures appearing in the vortex struc- 
ture, such as the pair potential, the internal magnetic 
field and the local density of states (LDOS). 

Low energy excitations, which dominantly govern the 
vortex structure, can be divided conceptually into three 
categories; (1) those from the continuum states associ- 
ated with the nodal structure in the d^2_y2-waYe pairing, 
(2) the core excitations from the bound states localized in 
a vortex core, and (3) the quasiparticle transfer between 
vortices, i.e., the vortex lattice effect. While marry cal- 
culations for the (ij,2„j,2-wave pairing such as Rcf. take 
into account only item (1), which is valid near the lower 
critical field Hd, items (2) and (3) are also indispensable 
when considering the vortex state in general. In order 
to help establish the general features of the mixed state 
in both the dx2-y2- and s-wave cases, one needs to cal- 
culate the vortex structure while taking into account all 
three factors on an equal footing. Through these efforts 
we may gain a more valid and vivid picture of the vortex 
for the whole region of Hd < H < Hc2 {Hc2 is the upper 



critical field). 

Detailed properties of the vortex structure, such as 
its field dependence, attract much attention even in the 
conventional s-wave superconductors. Several important 
means to probe the vortex structure are now available ex- 
perimentally in various superconductors, including high- 
Tc superconductors. Specific heat experimentstTEl in- 
form ija-pf low energy excitations. Muon spin resonance 
(/iSR)Lra and small angle neutron scattering (SANS) in- 
vestigate the internal field distribution of the-vortex 
structure. Scanning tunneling microscopy (STM)BtIj di- 
rectly observes the LDOS. These data are often analyzed 
within conventional phenomenological theories such as 
the Ginzburg-Landau (GL) theory or the London theory. 
The GL theory is, strictly pSoeaking, valid only near the 
transition temperature TcEj As for the London theory, 
which is applied near Hd^ the simple cutoff procedure 
for the core radius is too rough an approximatio n j j O ies- 
timate the exact contribution of the vortex core.tSEj In 
this sense it is necessary to develop a microscopic theory 
in order to correctly analyze valuable experimental data. 

In this paper, we investigate various aspects of the vor- 
tex lattice structure based on the quasi-classical Eilen- 
berger theory to clarify the difference in the vortex struc- 
ture between the dr^2_y2- and s-wave pairings. The case 
of a magnetic field applied along the c axis (or z axis) 
is considered. We follow the method of Rcfs. 16 and [T^, 
where the s-wave pairing case was investigated. This 
method can be applied in most regions of the mixed 
state, i.e., < T < Tc and Hd < H < Hc2- We cal- 
culate both the pair potential and the vector potential 
self-consistently in the vortex lattice case. Then we dis- 
cuss the spatial variation of the current and the internal 
field in detail. As for the single vortex case, the fourfold- 
symmetric vortex core structure of the d^2_y2--wave pair- 
ing was shown in Ref. |l^ by the quasi-classical calcula- 
tion. In the present paper, where we examine the vortex 
lattice case, we can see how this c?2,2_j,2-wave nature and 
the vortex lattice effect affect each other. And we also 
study the field dependence of the vortex structure. 

The rest of this paper is organized as follows. In Sec. 
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U, we describe the method of the quasi-classical calcula- 
tion. In Sec. in, we estimate the free energy for several 
configurations of the vortex lattice. There, the square 
lattice has lower free energy than the conventional trian- 
gular lattice at higher fields. In Sec. we study the 
field dependence of the vortex structure (the pair poten- 
tial, the current and the internal field) in the square lat- 
tice case. The form factor and the distribution function 
of the internal field are discussed in connection with the 
SANS and /iSR experiments. The vortex structure in the 
triangular lattice case is briefly discussed in Sec. We 
study the LDOS structure and the field dependence of 
the spatially averaged DOS in Sec. and the structure 
of th e ind uced order parameter of the other symmetry in 
Sec. Vn. A summary and discussions are given in Sec. 



VIII. 



II. QUASI-CLASSICAL EILENBERGER THEORY 

Our calculations are performed following the method 
of Refs. 16 and |l^. We consider the case of the clean limit 
and the cylindrical Fermi surface. These are appropriate 
to high-Tc superconductors and low-dimensional organic 
superconductors. 

First, we obtain the pair potential and the vector po- 
tential self-consistently by solving the Eilenberger equa- 
tion in the Matsubara frequency cun = (2n + 1)ttT. We 
consider the quasi-classical Green functions g{iujn,9,r), 
f{iujn, r) and f''{iujn, r), where r is the center of mass 
coordinate of the Cooper pair. The direction of the rela- 
tive momentum of the Cooper pair, k = k/|k|, is denoted 
by an angle 9 measuxed from the x axis. The Eilenberger 
equation is given byO 



{cj„ + %F • (- + ^ A(r)) }f{i^n, 0, r) 
= A(^,r)g(zo^„,0,r), 

{c.„-^VF-(|-^A(r))}/t(zc.„,0,r) 



/(zc^„,0,r)/t(jc^„,0,r)]i/2, 



(2.1) 

(2.2) 
(2.3) 



where Re(7(ia;„, 0, r) > and vp — wpk is the Fermi 
velocity. The vector potential is written as A(r) — x 
r -|- a(r) in the symmetric gauge, where H = (0, 0, H) 
is a uniform field and a(r) is related to the internal field 
h(r) = (0, 0, /i(r)), where h(r) = Vxa(r). As for the pair 
potential A(6', r) = A(r)0(6'), we set (t){6) = cos 2(61- 
6o) for the dx2_y2--wa,ve pairing and (j){0) — 1 for the s- 
wave pairing. Here, is the angle between the x axis and 
the a axis of the crystal coordinate. The self-consistent 
conditions for A(r) and a(r) are given as 



A{d, r) = iVo27rT ^ / —¥{9', 0)/(zc^„, 9', r) 



(2.4) 
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— -g[iujn,0,r), 
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(2.5) 



where A^o is the density of states at the Fermi sur- 
face, V{9',9) = V4){9')(t){9) is the pairing interaction, 
K = (7C(3)/72)i/2(Ao/rc)KBCS with Rieman's zeta func- 
tion C(3)- kbcs is the GL parameter in the BCS the- 
ory, and Ao is the uniform gap at T = 0. We set the 
energy cutoff ojc = 20Tc. In the following calculations, 
energies and lengths are measured in units of Aq and 
Co = i'f/Ao = TT^BCS (Cbcs is the BCS coherence length), 

respectively. 

By solviHig Eqs. ( |2.l| )-( ^^ by the so-called explosion 
methocOliZi under A(r) and a(r) of the vortex lattice 
case, we estimate the quasi-classical Green functions at 
40 X 40 discretized points in a unit cell of the vortex lat- 
tice. We obtain a new A(r) and a(r) from Eqs. (2.4) 
and (2.5), and use them at the next step in calculating of 
Eqs. ( |2.l[ )-( pT3| ). This iteration procedure is repeated un- 
til a sufficiently self-consistent solution is obtained. We 
use the material parameters appropriate to YBCO, i.e., 
$BCS = 16A and kbcs=100. There, i/c2=66.7 Tesla in the 
s-wave pairing and_93.2 Tesla in the (ia.2_j^2-wave pair- 
ing for r/Tc = 0.5.EJ To study the field dependence, the 
calculations are done for various fields at a fixed temper- 
ature T/Tc ~ 0.5. The spatial variation of the internal 
field and the current J(r) — (c/47r)V x h(r) is calculated 
from a(r). 

Next, we calculate the LDOS for energy E by 



N{E, r) = A^o 



—Re g[i^n 
At: 



E + nj,9,r). (2.6) 



To o btain g{iujn E + iri,9,r), we solve Eqs. (2.1)- 
( ^.3| ) for 1] — iE instead of ujn using the self-consistently 
obtained A(r) and a(r). The DOS is given by 



N{E)^{N{E,r)). 



(2.7) 



Here, (• • •) = /^^.^ ^^jj dr{- ■ ■)/S means the spatial aver- 
age, where S is the area of a unit 
Free energy F is calculated byli^ 



F 



, (Mr)^) (|A(r)p) 
(0o/eo')2 ^ NoVAl 
2^T ^ [^^ d9 



'0 a;„>0"'0 



(2.8) 



with 



/(^„,0,r) = A*(0,r)/ + A(0,r)/t 

-<'-4}{-+r-(l+>0}/ 

If i /V 27r 



A*(0,r)/ + A(0,r)/t 



1 + 5 



(2.9) 
(2.10) 
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where g, / and mean g{iu:n,0,r), fjiojn 9,r) and 



/^(iw„, 0, r), respectively. We use Eqs. ( ^.lD -( ^.3[) to 
obtain Eq. ( |2lo| ). 



III. VORTEX LATTICE CONFIGURATION 

We have to fix the shape and the orientation of the 
vortex lattice before calculation. The shape of the vortex 
lattice is characterized by the unit vectors ri = (a^;, 0, 0) 
and r2 = {Co-x, Q-j/, 0) where a^ciy = (fio/H and where the 
X axis is set in the ri direction. Since we consider the 
case |ri| — |r2|, we can write r2 — {ux cos9l, clx sinOL), 
where 6*^ is the angle between ri and r2. 

In this paper, we study the contributions of the d^^-y^- 
wave effect and the vortex lattice effect on the vortex 
structure. As the d^2_y2-wave nature reflects the node 
direction of the energy gap, its directional dependence is 
determined by the angle from the crystal coordinate. To 
investigate the dj.2_y2-wave nature, we compare the vor- 
tex structure along the a axis and b axis directions (we 
denote these directions as "the 0° directions" in the fol- 
lowing) and that along the 45° directions from the a axis 
(we denote these as "the 45° directions"). On the other 
hand, to study the vortex lattice effect, we compare the 
vortex structure along the nearest neighbor (NN) vor- 
tex direction and that along the next nearest neighbor 
(NNN) vortex direction. We denote the former direc- 
tion as the "NN direction", and the latter direction as 
the "NNN direction" . The orientation of the vortex lat- 
tice, i.e., the relative angle between the "0° directions" 
and the "NN direction" is important in considering the 
vortex lattice structure in the d^2_y2--w8ive pairing. The 
orientation is characterized by 9q. 

In an isotropic s-wave superconductor, the shape of 
the vortex lattice is a 60° triangular lattice. However, 
in a dx2-y2-wave superconductor, the shape may devi- 
ate from a triangular lattice. The deviation comes from 
the fourfold-symmetric vortex structure around the vor- 
tex core, which reflects the d^2_y2-wave symmetry of the 

pair potential in the k space (i.e., — ky). To deter- 
mine the vortex lattice conflguration of the equilibrium 
state, we compare F of Eq. ( ^.8|) for the various vortex 
lattice configurations, and find the state having mini- 
mmn F. As our quasi-classical calculation needs a lot of 
time to compute, we cannot check all possible configu- 
rations of shape and orientation. So, we only compare 
the square lattice case (0l = 90°) and the triangular lat- 
tice case {Ol = 60°) in the dr^2_y2-waMe pairing. There, 
we consider the two orientation cases 0^ — 45° and 0°. 
For 6*0 = 45° (6*0 =0°), one of the NN vortices is lo- 
cated in the 45° (0°) direction. The field dependence of 
F is shown in Fig. |l|. There, we present the difference 
F — i^sq45, where i^sq45 is F for the case of the square 
lattice with 6q = 45°. The square lattice with 9q — 45° 
has the smallest F at a higher field H/Hc2 > 0.15. At 
a lower field, the triangular lattice case has the smallest 



F. So, the vortex lattice configuration changes from the 
triangular lattice to the square lattice with 9q = 45° as 
the field increases. We note that there may be a sta- 
ble oblique lattice between the triangular lattice and the 
square lattice. An estimate including the oblique lattice 
by the quasi-classical theory belongs to a future work. As 
the square lattice with 0o = 0° has a higher F, it is an 
unstable state. As for the orientation of the triangular 
lattice, the case 9q = 45° has a smaller F than the case 
6*0 = 0°, but the difference is very small. 

Some theoretical calculaticps i-epnsidered the vortex 
lattice confiauxation near Hc2 ,EirE3 or based on the Lon- 
don theory.EirEa These theories suggested that the vortex 
lattice configuration gradually changes from the triangu- 
lar lattice to the square lattice with increasing field, and 
the oblique lattice is realized at a lower field before set- 
tling into the square lattice. The square lattice config- 
uration with Of) = 45° has lower free energy than that 
of the conventional 60° triangular lattice over a wide re- 
gion of higher fields and lower temperatures in the mixed 
state for the dx2-y2-^ave pairing. This is consistent 
with the estimate of Fig. [T]. An STM experiment on 
YBa2Cu307_5 (YBCO) also suggested the squaKe,vortex 
lattice configuration in the c?2.2_j,2-wave pairing .E2J 



IV. VORTEX STRUCTURE IN THE SQUARE 
LATTICE CASE 

As noted in the previous section, the stable configu- 
ration of the vortex lattice at higher field is the square 
lattice with 6*0 = 45°. We fix the shape of the vortex 
lattice as the square lattice in order to study the field 
dependence of the vortex structure. To clarify the roles 
of the (i2,2_y2-wave effect and the vortex lattice effect, 
we present the results of the unstable orientation case 

= 0° in addition to the stable orientation case 6*0 = 45° 
in the (i^2_y2-wave pairing. We denote the former as "the 
c?0° case" , and the latter as "the (i45° case" . The s-wave 
pairing case is also calculated for the square lattice to 
clarify the effect of low lying excitations associated with 
the gap anisotropy. The comparison among these cases 
(s-wave case, dAb° case and dQ° case) helps us to under- 
stand the contributions of the (ia;2_y2-wave nature and 
the vortex lattice effect. 



A. Amplitude of the pair potential 

In Fig. ^, we show the spatial variation of the vortex 
structure at low field H/Hc2 = 0.021 for the dAh° case. 
We show their contour lines and their profile. The profile 
is presented along the lines OA, OB, OC and AC shown 
in Fig. |. The line OA (OC) is the radius in the NN 
(NNN) direction, and the line AC is along the boundary 
of the Wigner-Seitz cell of the vortex lattice. 
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The amplitude |A(r)| is presented in Figs. || (a) and 

(d) . At low field, the vortex core region occupies a small 
part of a unit cell of the vortex lattice. There, |A(r)| 
approaches a constant value (0.95Ao for T = 0.5Tc) at a 
few ^0 away from the vortex center, as presented in Fig. 
H (d). The d^, 2 _ J, 2 -wave nature appears in the shape of 
the contour lines of |A(r)|. In the s-wave pairing, the 
contour lines are circular around each vortex center at 
low field. In the d3,2_j,2-wave pairing, as shown in Fig. |^ 
(a), the contour lines become a fourfold-symmetric shape, 
especially for O.SAq < |A(r)| < 0.9Ao. There |A(r)| is 
suppressed around the vortex core in the 0° directions 
compared with the 45° directions. This structure is con- 
sistent with that obtained by the single vortex calcula- 
tion.El 

Figure ^ shows |A(r)| at a higher field H/Hc2 = 0.54 
for three cases. At high field, the vortex core region occu- 
pies a large part of a unit cell. So, |A(r)| does not reduce 
to a uniform value even in the boundary region, as shown 
in Figs. ^ (d)-(f). The influence of the vortex lattice ef- 
fect on |A(r)| can be seen in Figs. || (a) and (d) for the 
s-wave pairing case. While the vortex core has a circu- 
lar symmetric structure in the inner region of the core 
(|A(r)| < 0.5Ao), |A(r)| shows directional dependence 
in the outer region, reflecting the NN vortices. There, 
|A(r)| is suppressed in the NN directions (line OA) com- 
pared with the NNN directions (line OC). In the dj;2_y2- 
wave pairing, the contribution of the fourfold-symmetric 
core structure appears in addition to the vortex lattice 
effect, depending on the vortex lattice's configuration, 
and especially on its orientation. We demonstrate this 
by showing two orientation cases d45° [Figs. I (b) and 

(e) ] and dO° [Figs. || (c) and (f)]. In both cases, the in- 
ner region of the vortex core shows a fourfold-symmetric 
structure, seen also in the low field case shown in Fig. ^ 
(a). However, the difference between the two orientation 
cases appears in the outer region of the vortex core. Due 
to the (ia;2_j,2-wave nature, |A(r)| is suppressed in the 
0° directions. Because of the vortex lattice effect, |A(r)| 
is suppressed in the NN directions. When these two di- 
rections disagree as in Fig. ^ (b), these two suppres- 
sions cancel each other. The difference between |A(r)| 
in the NN direction (line OA) and |A(r)| in the NNN 
direction (line OC) is smeared. Along the boundary line 
AC, |A(r)| is almost constant. On the other hand, when 
the 0° direction and the NN direction agree as in Fig. 
^ (c), the suppressions due to the vortex lattice effect 
and the (i^2_y2-wave effect enhance each other. |A(r)| 
is largely suppressed in the 0° directions (i.e., the line 
OA in the NN direction) compared to |A(r)| in the 45° 
direction (i.e., the line OC in the NNN direction). The 
result shown in Fig. |l], which shows that the (i45° case is 
stable at higher field and the dO° case is unstable, may 
be related to the difference in the vortex structure shown 
in Figs. I (b) and (c). 

Next, we study the field dependence of |A(r)|. We de- 
fine Ann (Annn) as the maximum of |A(r)| along the 
line OA (OC) in the NN (NNN) direction. The field de- 



pendence of Ann and Annn is presented in Fig. || for the 
s-wave (a) and the dj.2_y2-wave (b) pairings. They de- 
crease to as Hc2 is approached. The difference between 
Ann and Annn, which refiects the vortex lattice effect, 
appears when H/Hc2 > 0.1. The difference is small for 
the (i45° case, but large for the dQ° discussed 
above. 

The field dependence of the core radius is shown in 
Fig. ^. The radius is defined from the initial slope of 
the pair potential by setting as |A(r)| = AnnJ'/Ci at the 
vortex center. As H increases, the radius decreases 
similarly for both the s-wave and the d45° cases. In the 
unstable dO° case, the decrease of is saturated at a 
higher field compared to the other cases. 

B. Internal current distribution 

The spatial variation of the current |J(r)| at low field 
H/Hc2 = 0.021 is shown in Figs. | (b) and (e) for the 
dAb° case. There, |J(r)| has four peaks around each vor- 
tex core in the 45° directions at r ~ 0.5^o- The four 
peaks are dx^-y^-wtcve natuce, and are consistent with 
the single vortex calculation.ES In the s-wave pairing at 
low field, |J(r)| distributes circularly around each vor- 
tex without four peaks. Outside the vortex core region 
(r > ^o): there is no clear difference between the s- and 
(i2.2_j^2-wave pairings in the current distribution. 

Figure]^ shows | J(r)| at high field H/Hc2 = 0.54, which 
corresponds to the case of Fig. ^. In the s-wave pairing 
case of Fig. ^ (a), there are four peaks around each vortex 
core in the NN vortex directions. These peaks are due 
to the vortex lattice effect and enhance the four peaks of 
the dj.2_y2-wdMe nature in the 45° directions for the dAb° 
case [Fig. ^ (b)], but smear them for the dQ° case [Fig. 
0(c)]. 

We can also define the core radius by the maximum 
current. The core radius ^2 is defined as the radius where 
|J(r)| is at a maximum in the NN vortex direction. The 
field dependence of ^2 is presented in Fig. || (b). The 
radius ^2 decreases similarly to ^1 as shown in Fig. ^ (a). 
From both definitions of ^1 and ^2 (i-e., the initial slope 
of |A(r)| and the maximum current), we obtain a similar 
dependence on H. The shrinkage of the corepaidius was 
also reported in STM and ^SR experiments. li^u There, 
the results of the field dependence were analyzed in the 
dirty limit theory by the Wigner-Seitz method (i.e., a 
circular unit cell is used instead of the periodic boundary 
condition of the vortex lattice) . In our study, we obtain 
the field dependence in the clean limit, where we treat 
the periodic boundary condition exactly. 

The radius is also presented for the triangular vortex 
lattice of the s-wave pairing[Fig. ||. The field depen- 
dence is almost the same as in the square lattice case of 
the s-wave pairing. We confirm that the core radius of 
the d^2_y2-w&ve pairing shows a similar decrease to the 
triangular vortex lattice case[Fig. ^. When we calculate 
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the radius for a smaller kbcs ('^bcs = 7) in the s-wave 
pairing to study the case of NbSe2 , a similar decrease is 
obtained. So, the decrease of the core radius occurs over 
a wide range of kbcs- 



C. Internal field distribution 

Reflecting the four peaks of |J(r)| in the 45° directions, 
the internal field distribution has a fourfold-symmetric 
structure around each vcptex core, as presented by the 
single vortex calculation.ta h{r) extends toward the 0° 
directions in which the screening current is weak. How- 
ever, this fourfold-symmetric structure is restricted to 
the vortex core region, since the (ia;2_y2-wave nature is 
compensated by J(r) outside the core. Outside the core 
region, /i(r) immediately reduces to a circular structace- 
This compensation is seen in the single vortex caseJl3'EZI 
but it also appears in the vortex lattice case at low field. 
In Fig. H (c) and (f), we show the spatial variation of 
h{r) at low field H/Hc2 = 0.021 for the rf45° case. In the 
da;2_j,2-wave pairing, the contour lines of h{r) in the vor- 
tex core region are distorted from the circle of the s-wave 
pairing case into a fourfold symmetry. The distortion can 
be seen within the small core region as shown in Fig. |^ 
(c). Outside the core region, h{r) recovers the circular 
contour lines seen in the s-wave case. So the dx2_y2-weLYe 
nature scarcely affects the internal field distribution at 
low field. In the s-wave pairingcase, h{r) has the same 
structure as that shown in Fig. g (c) except for the small 
region of the vortex core. 

Figure || shows h{r) at high field iJ/_ffc2=0.54, which 
corresponds to the case shown in Fig. |^. In the s-wave 
pairing, h{r) is lowest at the boundary of the NNN di- 
rection (i.e., the point C in Fig. ^). There is a saddle 
point in h{r) at the boundary of each NN direction (i.e., 
the point A in Fig. ^). The fourfold symmetry around 
the vortex core due to the d^2_,,2-wave nature is seen in 
Figs, g (b) and (c), where /i(r) is enhanced in the 0° 
directions, and suppressed in the 45° directions. In the 
d45° case in Fig. ^ (b), the c?a.2_y2-wave nature enhances 
h{v) at the minimum point in the NNN direction, and 
suppresses it at the saddle point in the NN direction. So, 
variation of h{v) along the boundary line AC is small. 
On the other hand, for the dO° case in Fig. ^ (c), h{v) 
is enhanced at the saddle point and suppressed at the 
minimum point due to the (ia;2_y2-wave nature. So, the 
difference between the field at the saddle point A and 
the field at the minimum point C is large. Since the /i(r) 
distribution is observed in SANS and ^SR experiments, 
we discuss the related quantities in detail. 

First, we investigate the form factor of h{r). It is mea- 
sured by the SANS experiment. The form factor hm,n 
(to and n are integers) is the Fourier component of h{r) 
defined as 



where the reciprocal lattice vector gm.n = —nk.i -\- mk2, 
where ki — 2ht2 x t^/S, and k2 = 27rr3 x Yi/S with 
ra = (0, 0, 1) and S" = ri • (r2 x r^) ■ The field dependence 
of the dominant form factor hi o is plotted in Fig. ^ 
for the s- and dj,2_j,2-wave pairings. The two pairing 
cases show a similar dependence on the applied field. So, 
it is difficult to detect the c?2,2_j,2-wave nature in /ii,o- 
However, we can detect it in the higher order components 
such as /i2,o and ft-i^. The field dependence of /i2,o//ii.Oi 
/ii,i/^i,Oj ^3,o/^i,o and ft.2,i//ii,o is plotted in Fig. |l^. 
The effect of the d3.2_j,2-wave nature appears depending 
on the orientation of the vortex lattice. As shown in 
Fig. [l^ (a), with increasing field, ft.2,0/^1,0 of the s-wave 
pairing decreases and becomes negative. In the stable 
d45° case, /12, 0/^1,0 remains positive. However, in the 
unstable dO° case, /i2. 0/^1.0 decreases more rapidly than 
that of the s-wave pairing. As shown in Fig. 10 (b), 
^1,1/^1,0 remains almost constant at high field in the 
s-wave pairing. In the d45° case, ft.i,i/ft.i,o increases as 
Hc2 is approached. But, in the c?0° case, it decreases. For 
6'o = 0° and for = 45°, the d^2 _y-i-^scve effect shifts 
the form factors in the opposite direction. The d^^^y^- 
wave nature also appears in /13, 0/^1,0 [Fig- (c)] and 
h2,i/hi,Q [Fig. ^ (d)]. We expect that these d^-^-y^- 
wave natures can be detected by the SANS experiment. 
The form factors are also presented for the triangular 
vortex lattice of the s-wave pairing in Figs. || and |l^. 
Their field dependence is qualitatively the same as in the 
square lattice case of the s-wave pairing. 

Next, we investigate the magnetic field distribution 
function defined as 



P{h) = {5{h-h{v))). 



(4.2) 



/i(r) ^ hm,n exp(ig„i,„ • r) 



(4.1) 



It corresponds to the resonance line shape in /iSR or 
NMR experiments. In Fig. ^ we plot P{h) for 
i//7?c2=0.021(a), 0.15(b) and 0.54(c) in the s-wave pair- 
ing case and in the two orientation cases of the d^^-y^- 
wave pairing. We denote the maximum (minimum) of 
h{v) as /imax (^min)- Reflecting the saddle point in /i(r), 
P{h) has a logarithmic singularity at hg. In the low field 
case shown in Fig. ^ (a), P{h) is almost the same shape 
for both the s- and o?2,2_y2-wave pairings. For compari- 
son, we also present the line shape of the triangular lat- 
tice case in Fig. |ll| (a), which also does not show any 
difference between the s- and dj.2_j,2-wave pairings. As 
shown in Fig. ^ (c), the d^2 _y2-^ffsve nature scarcely ap- 
pears in /i(r) at a low field. With increasing field, the 
dx2 -y2-wa^re nature gradually appears in P{h) as shown 
in Figs. |l^ (b) and (c). For the stable c?45° case, /imin in- 
creases and hg decreases cornpared with the s-wave pair- 
ing case, as shown in Fig. || (b). For the unstable (iO° 
case, /imin dccrcascs and hg increases. These changes of 
/iinin and hg can be seen in the resonance line shape in 
Figs. ^ (b) and (c). We note that the line shape of the 
d45° case in Fig . [ll| (c) resembles that of the triangular 
lattice (see Fig. |ll| (a)). So, in the d3;2_j,2-wave pairing, 
it is risky to determine the vortex lattice shape by the 
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conventional method based on the distance between hs 
and /imin, i.e., (/is - /imin)/(/imax " /imin)- The field de- 
pendence of /imax, hs and /imin IS presented in Fig. 
There, /imax shows a similar dependence for both the s- 
and dj,2_^2-wave pairings. We also calculate the field de- 



in the triangular vortex lattice for the 
j,2-wave pairings. Their dependence is 



pendence of ft-max 
s-wave and d^^ 
almost the same as that of the square lattice case. So, 
/imax is not much affected by the vortex lattice configu- 
ration and the pairing symmetry. The dj:2_y2-wave na- 
ture affects /imin and hg. For the stable (i45° case, /imin 
approaches hg with increasing field, and /imin ^ /is for 
H/Hc2 > 0.5. We hope that this characteristic of the 
resonance line shape P{h) will be examined by the /iSR 
or NMR experiments. 

To help in analyzing the shape of the distribution 
hmction P{h), we plot the field dependence of the vari- 
ance a = (/i(r)^)^/^ and of the skewness parameter 
a — (/i(r)'^)^/^/cr in Fig. l^j-JCijejsa quantities are ana- 
lyzed by the /xSR experiment £j'E2ti23 As for the variance, 
both the s- and dx2_y2-wave pairings show a similar de- 
pendence on H. But, the difference between the s- and 
dx2-y2-weLve pairings clearly appears in the skewness pa- 
rameter. 



V. VORTEX STRUCTURE IN THE 
TRIANGULAR LATTICE CASE 

To study the vortex structure at low field and the effect 
of the vortex lattice shape, we also consider the triangular 
vortex lattice case. In Fig. Q we show the vortex struc- 
ture at H/Hc2 = 0.15 for the dj.2_j^2-wave pairing. The 
contour lines of |A(r)|, |J(r)| and /i(r) are presented for 
the two orientation cases 9q = 45° and 0°. The contour 
lines of |A(r)| show a fourfold-symmetric shape, |J(r)| 
has four peaks in the 45° directions, and /i(r) extends 
in the 0° direction around each vortex core. These fea- 
tures show the dx2-y2-waNe nature presented in the pre- 
vious section. At a low field such as H/Hc2 = 0.021, 
the fourfold-symmetric structure of the (i^2_j^2-wave na- 
ture is restricted to the small region of the vortex core 
as in the case shown in Fig. |[ These core structures 
are not affected by the shape and the orientation of the 
vortex lattice at low field. But, with increasing field, 
the fourfold-symmetric core structure gradually feels the 
effect of the hexagonal Wigner-Seitz cell of the vortex 
lattice. By the effect of the fourfold-symmetric struc- 
ture around the core region, the six NN directions of the 
triangular lattice are not equivalent in the dx2^y2-wave 
pairing. In Fig. |lj (a) for Oq = 0°, the suppression 
of |A(r)| in the NN directions due to the vortex lattice 
effect is enhanced in the a axis direction by the d^2_y2- 
wave nature. So, the structure of the two NN directions 
along the a axis is not equivalent to that of the other four 
NN directions. The suppression due to the d3;2_j,2-wave 
nature in the b axis direction is weakened by the vortex 
lattice effect because it is in the NNN direction. The 



vortex structure now shows twofold symmetry, which be- 
comes more prominent at higher field. In Fig. (d) for 
6*0 = 45°, the suppression in the NN direction is weak- 
ened in the 45° direction, compared with the other four 
NN directions. 

The structure of these two NN directions can also be 
seen in the current |J(r)|jFigs. |l^ (b) and (e)] and the 
internal field h{r) [Figs. |l| (c) and (f)]. In Fig. |l| (c) 
for 9o = 0°, h{r) at the two saddle points in the a axis 
NN direction is larger than that of the saddle points in 
the other four NN directions. In Fig. |l^ (f) for Oq = 45°, 
/i(r) at the two saddle points in the 45° direction becomes 
small compared to the other four NN directions. The 
resonance line shape P{h) in the triangular lattice case 
is presented in Fig. where the line shape around the 
peak is focused. Refiecting the two saddle points of /i(r) 
in Figs. ^ (c) and (f), the line shape of P{h) has two 
logarithmic peaks in the da;2_j^2-wave pairing. The two 
peak structure was suggested for the oblique lattice state 
between the square lattice and triangulaij- 
on the GL theory or the London theory.ol 
occurs even in the triangular lattice for the dj.2_y2-wa,ve 
pairing by the influence of the fourfold-symmetric vortex 
core structure. 

Relaxing the difference between these two NN direc- 
tions, the shape of the vortex lattice may deviate from 
the 60° triangular lattice. 



based 
But, it 



For 6*0 = 0° in Fig. |14| (a), 
the intervortex distance in the a axis NN directions may 
increase compared with the other NN directions. For 
6*0 = 0° in Fig. |lj (d), the intervortex distance in the 
45° directions may decrease compared with the other NN 
directions. This deformation may relax the difference be- 
tween the two NN directions. This is the origin of the 
distortion from the conventional 60° triangular lattice. 
The oblique lattice is realized in the intermediate field 
region between the square lattice and the 60° triangular 
lattice. In the oblique lattice, the double-peaked struc- 
ture of P{h) shown in Fig. 



15 may be smeared, since the 



difference of the two saddle points is weakened by the 
distortion from the 60° triangular lattice. 

We also calculate the field dependence of the form fac- 
tor hm.n in the triangular lattice case. It shows a similar 
dependence to the square lattice case presented in Figs. 
^ and |l^. The dx2-y2-wa,Ye nature does not much affect 
hi Q, but appears in the higher order factors such as /ii^i 
and /i2 • 



VI. DENSITY OF STATES 

The LDOS N{E,r) is calculated by Eq. (|2j) in the 
square vortex lattice case for the s-wave and the d45° 
cases. Typically, we choose rj = 0.03. First, we dis- 
cuss the LDOS for the s-wave pairing. In the calcula- 
tion for a single vortex, the LDOS for low energy exci- 
tations is mainly distributed in a circle around the vor- 
tex center Ej'E^ The radius of the circle increases if E is 
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raised. The LDOS of the vortex lattice case is presented 
in Fig. ^ for the s-wave pairing at H/Hc2 = 0.15. The 
area presented is the Wigner-Seitz cell of the square vor- 
tex lattice, i.e., within the dotted line shown in Fig. 
The LDOS for low energy excitations is distributed in a 
circle as in the single vortex case. Further, by the vortex 
lattice effect (i.e., the quasiparticle transfer between vor- 
tices), N{E = 0, r) is suppressed on the hue connecting 
two NN vortex centers, as shown in Fig. ^ (a). A small 
suppression also occurs along the line between two NNN 
vortex centers. For a finite E below Aq, the suppression 
occurs in the NN directions along the tangent lines of 
the circle of the dominant LDOS distribution, as shown 
in Fig. |l^ (b) and (c). The suppression is scarcely seen 
at low field {H < 0.1Hc2)- There, the LDOS distribution 
shows an almost circular symmetry around each vortex. 
At high field, however, the suppression becomes signif- 
icant and the LDOS distribution deviates from the cir- 
cular symmetry and shows fourfold- (sixfold-)symmetric 
structure in the square (triajagular) lattice case, reflect- 
ing the vortex lattice effect.Ej At high field, the LDOS 
has a finite distribution all over the unit cell in addition 
to the ridge structure, while at low field it is distributed 
in a restricted region along the circular ridge. 

Next, we discuss the d^2_y2-wave pairing case. The 
single vortex calculation shows that the LDOS for low 
energy excitations around a vortex has prominent tails 
extending in the 45° directions [see Fig. 8 (d) of Ref. 

reflecting the node structure of the energy gap in 
the d2;2_y2-wave pairing. Since the low energy quasipar- 
ticle distribution around a vortex extends toward infi- 
nite points in the 45° directions, this low eaergy state is 
not a bound state in the exact meaningEil The LDOS 
of the vortex lattice case is presented in Fig. |l^ at 
H/Hc2 = 0.15 for the (i45° case. In this configuration, 
the 45° directions are the NN vortex directions, which 
are on the x axis and y axis in the figure. In N{E = 0, r) 
in Fig. ^ (a), the tail structure in the 45° directions in 
the single vortex case is modified due to the vortex lattice 
effect. Each tail of the zero energy LDOS in the 45° di- 
rections splits into two ridges because of the suppression 
by the vortex lattice effect in the NN vortex directions. 
So, the tails extend in rather different directions from the 
45° directions. In the (ia;2_j,2-wave pairing, the suppres- 
sion of the vortex lattice effect occurs from a low field 
such as H/Hc2 = 0.02. It means that the quasiparticle 
transfer between vortices is large due to the tail struc- 
ture of the LDOS in the d^2_y2-wa,Ye pairing, compared 
with the s-wave pairing case. With increasing field, the 
vortex lattice effect increases and the splitting of the tail 
structure in the 45° directions becomes pronounced. 

For a finite E below the energy gap a/2Ao, N{E,r) 
shows a -shape ridge structure in the (i^2_j^2-wave 
pairing, as shown in Fig. |l^ (b) and (c). This is copsife 
tent with the calculation for the single vortex caseJl3'E2l 
The ridge structure is not much affected by the vortex 
lattice effect except when E ^ 0. The reason is as fol- 
lows. The quasiparticle flow is along the ridges on the 



four open trajectories.tS When the flow direction 9 of 
the quasiparticles on the trajectories approaches the 45° 
directions, it leaves the vortex center. There, the quasi- 
particles with a finite E become a scattering state, since 
y/2Ao\(j){e)\ < Eis satisfied. When the quasiparticle flow 
is toward the NN vortex directions, the vortex lattice ef- 
fect occurs. However, since the NN directions are now 
the 45° directions, the quasiparticles flowing toward the 
NN vortex directions on the trajectories are in a scatter- 
ing state. A scattering state is widely distributed and 
docs not contribute to the ridge structure. 

As a result, the local spectrum N{E, r) is given as fol- 
lows in the vortex lattice state. There are some peaks 
below the energy gap as in thp-single vortex case, reflect- 
ing the trajectory structure.E3 Extra small peaks may 
also appear by the trajectories extending from neighbor- 
ing vortices. And the low energy state for i? '--^ is 
suppressed due to the vortex lattice effect on the line 
connecting two NN (or NNN) vortex centers. 

We have considered the case of the square lattice with 
6*0 = 45°. From that information, we can estimate the 
spatial variation of N{E = 0, r) in the other vortex lat- 
tice configuration. First, we consider the contribution 
on N{E — 0, r) from the vortex core and the four tails 
extending in the 45° directions from each vortex center. 
The suppression along the line connecting NN vortices is 
also taken into account. There is also a small suppres- 
sion in the NNN vortex directions. The suppression is 
larger at higher field. Then, we obtain the N{E = 0, r) 
structure in the other vortex lattice configuration. 

The relation N{E = 0) cx is sometimes considered 
for the c?a;2_y2-wave pairin|g. The relation is based on the 
single vortex calculation,EI which yields a structure with 
four tails extending in the 45° directions. These result 
from the node structure of the energy gap. However, 
since the vortex lattice effect smears the tail structure, 
the relation N{E = 0) (x ^/H may be affected. We ex- 
amined it in Ref. ^ and obtained N{E = 0) cx H^ *^. 
The exponent becomes smaller than 1/2. The relation 
N(E = 0) oc H is sometimes considered for the s-wave 
pairing. In this case, N{E = 0) is proportional to the 
vortex density since the low energy excitation consists of 
a bound state around each vortex core. We also examined 
this case in Ref. | and obtained N{E = 0) cx H°-^'^. The 
deviation from the iJ-linear relation results from the field 
dependence of the vortex core radius which we presented 
in Fig. ^. 



The spatially averaged DOS in Eq. (2.7) is presented 
in Fig. ^ for the s-wave (a) and (2x2 _j^2 -wave (b) pair- 
ings. We note that the energy gap in the (i^2_j^2-wave 
case IS V2Ao in Fig. |l|(b). The energy gap structure at 
zero field is gradually buried by the low energy excitation 
of the vortex as the field increases. While the local spec- 
trum N{E, r) has several peaks depending on the pairing 
symmetry and the position,E3 there is no peak in N{E) 
below the energy gap. The peak structure is smeared 
when the spatial average is taken. 
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VII. INDUCED ORDER PARAMETER OF 
OTHER SYMMETRY 

So far, we have considered the case where we can ne- 
glect the induced order parameter of other symmetry, i.e., 
the pure d^2_j^2-wave case. When there is s-wave {dxy- 
wave) interaction in addition to the dominant (1^2 _y2- 
wave pairing interaction, the s-wave (dajy-wave) compo- 
nent of the pair potential is induced at the place where 
the dominant dj.2_j,2-wave pair potential spatially varies, 
such as a vortex or an interface. The induced s-wave 
component was mainly studied by the two-component 
GL theory, where the coupling of the s-wave and the 
dj.2_y2-wa,Ye components in the gradient term plays an 
important rolcEirEj The induced d^y-wave component 
can be discussed by extending the two component GL. 
theory by including the higher order gradient coupling.E^ 
The quasi-classical calculation for the induced s-wave and 

4y-wav&xpmponents was performed for the single vor- the amphtudes |As(r)| and \Axy{r)\ in Fig. ^ Around 
tex case.E3 The induced component around the vortex 
was also reported in a studji-of jtiie Bogoliubov-de Gennes 
theory for several models. oeZtEj There, the amplitude of 
the induced component depends on the interaction pa- 
rameter and the carrier density. 

In this section, we study the spatial structure of the 



those found in a single vortex.Es With increasing field, 
since the size of the unit cell decreases, the structures 
of |As(r)| and |Aa;y(r)| gradually change to those shown 
in Figs. |9| (c) and (d), respectively. These structures 
have extra zero points on the boundary in addition to 
the vortex center. 

To understand the structures of As(r) and Axy{r), 
we present the structure of the phase schematically 
for low and high field in Fig. There, the posi- 

tion of the singularity and its winding number are pre- 
sented for the relative phase to the dominant Aj,2_j^2(r), 
i.e., arg(As(r)/A^2_j,2(r)) and arg(Aa;.y(r)/A^2_y2(r)). 
Since Ax2_y2(r) has the singularity of the winding num- 
ber at the vortex center, the winding number of 
arg(As(r)) itself is —1 at the vortex center [the case 
shown in Fig. ^ (a) and (c)]. And that of arg(Aa;j,(r)) 
is -3 at low field [Fig. |^ (b)] and -1-1 at high field [Fig. 
pO| (d)]. The singularity corresponds to the zero point of 



induced s- and dxy-'wa.ve components in the vortex lat- 
tice of c?2;2_y2-wave superconductors, based on the quasi- 
classical theory. These induced components can be cal- 
culated using Eq. (2.4), by setting 



A(6»,r) ^ A,{r)(l3,{9) + Ax2^y2{r)(t)x2^y2{9) 

+Axyir)(bxy{0), (7.1) 

V{e\e) = VsM0')M0)+Vx--y2<j)x--y-{0>.-^y^{O) 

+ %y(b.y{0')(b.y{O) (7.2) 

where <j)s{d) = 1, cpx^-y^id) — V2cos2(6' — Oq) and 
4'xy{()) = \/2sm2{6 — ^o). Note that the induced com- 
ponents As(r) and A^yir) are, respectively, proportional 
to Vs and Vxy Here, we consider the case [t^j, |t4y| <C 
\Vx2^y2\. In this case, the contributions of As(r) and 
Axy{r) are small and can be neglected compared to the 
domi nant A^2_ „2(r) contribution in the calculation of 
Eqs. (^-(pl). 

Figure |19| shows the spatial variation of the ampli- 
tudes |As(iF)| and |A2,j^(r)| in the dA5° case at low field 
H/Hc2 = 0.021 and high field H/Hc2 = 0.54. The pre- 
sented area is the Wigner-Scitz cell of the square vortex 
lattice, i.e., within the dotted line in Fig. H At low field, 
the contour line of |As(r)| is a four-lobed shape, as shown 
in Fig. |l^ (a). There, |As(r)| has four peaks around the 
vortex core in the 45° directions (which correspond to 
the X axis and y axis in this figure). The contour line of 
|Aa;jy(r)| is eight-lobed shape shape, as shown in Fig. |l^ 
(b). There are eight peaks of |A2:j,(r)| around the vor- 
tex core in the 0° and 45° directions. These structures 
of |As(r)| and |A2;y(r)| at low field are consistent with 



the singularity of the winding number n, the amplitude 
increases as r'"' . 

At low field, as shown in Fig. |2^ (a), the relative phase 
of the s-wave component has a —2 singularity at the vor- 
tex center and four -|-1 singularities in the 0°,--directions. 
This is consistent with the single vortex case.tJ Since the 
winding number of the relative phase should total zero 
for each unit cell to satisfy the periodic boundary con- 
dition, there is a —2 singularity at each corner of the 
Wigner-Seitz cell. With increasing field, the -1-1 singu- 
larity approaches the boundary and the four -f 1 singu- 
larities combine with the —2 singularities at the corners 
producing +2 singularities, as shown in Fig. |2^ (c). In 
the triangular vortex lattice case, each +2 (—2) singu- 
larity at the corner of the boundary splits into two -|-1 
(—1) singularities which are located at the corner of the 
hexagonal Wigner-Seitz cell. 

On the other hand, arg(A2;j,(r)/A3,2_y2 (r)) has a —4 
singularity at the vortex center and eight -1-1 singulari- 
ties around the vortex coie^ low field. This is consistent 
with the single vortex case.L^I Further, to satisfy the peri- 
odic boundary condition, there are eight —1 singularities 
on the boundary, as shown in Fig. ^ (b). As the field 
increases, the -1-1 singularities approach the —1 singular- 
ities on the boundary, becoming -1-1 singularities on the 
boundary, as shown in Fig. ^ (d). The —4 singularity at 
the vortex center splits into four —1 singularities around 
the vortex core. 

We also consider the d0° case, where the 0° direction 
is along the x and y axes. The four-peaked structure of 
|As(r)| at low field rotates 45° from its position in the 
(i45° case [Fig. |l^ (a)]. The low field structure is fixed to 
the crystal coordinate. However, as the field increases, 
the structure of As{r) reduces to that shown in Fig. |l^ 
(c) . The high field structure is independent of the crystal 
coordinate and is fixed to the shape of the vortex lattice. 
On the other hand, the structure of Axy{r) is the same 
as in the (i45° case regardless of the field. 

Even in an s-wave superconductor, if the extra dx2-y2- 
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wave (dxy-wave) component exists in the interaction, 
A3;2_y2(r) {Axy{r)) is induced around the vortex. We 
also calculate its structure in the vortex configuration 
shown in Fis;. I (a). The induced A^2_y2{r) {A^y{r)) 
has the same structure as As(r) in the (i45° {dO°) case. 

The induced component may be di^ussed— qualita- 
tively by the two component GL theoryO'EllE^Ga There, 
the vortex structure was studied in the single vortex 
case near Hd and in the vortex lattice state near 
Hc2- The results agree with ours. Near Hc2, the vor- 
tex structure is described by the magnetic Bloch state 
functiop ?/;„(r|ro), where n represents the n-th Landau 
leveLu As for the dominant dx2_y2-wa,ve pair poten- 
tial, A^2_j^2(r) ~ -^olrlro) in the leading order. For 
the induced s-wave component, As(r) ~ '02(r|ro) since 
it is induced by the gradient coupling {d^ — dy)A^2_y2 
in the GL equation. For the induced da;j,-wave compo- 
nent, Axy{r) ^ ip4{r\ro) since it is induced by the gra- 
dient coupling dxdy{dl — dy)Ax2_y2 in the extended GL 
equation.E3 In fact, at high field, As(r) and A^y(r) as 
obtained by our calculation have the same structure as 
'(/'2(r|ro) and '(/'4(r|ro), respectively. 



VIII. SUMMARY AND DISCUSSIONS 

We have studied the vortex structure and its field de- 
pendence in the dj.2_j^2-wave and s-wave superconductors 
by the self-consistent calculation of the quasi-classical 
Eilenberger theory for the vortex lattice case. At a 
higher field, the square vortex lattice configuration with 

= 45° (i.e., the NN vortex is located in a 45° direction 
from the a axis) has lower free energy than the conven- 
tional 60° triangular lattice for the (i^2_j,2-wave pairing. 
Then, we investigated the field dependence of the vor- 
tex structure in the square lattice configuration. Due to 
the (i2,2_j^2-wave nature, the contour lines of the ampli- 
tude |A(r)| and the internal field /i(r) show a fourfold- 
symmetric shape around the vortex core. The fourfold- 
symmetric structure is prominent in the vortex core re- 
gion. At low field, the (ia;2_y2-wave nature scarcely affects 
the structure of the outer region of the core. With in- 
creasing field, the core region occupies a larger area of 
the unit cell in the vortex lattice. Then, the d^2_y2-wave 
nature gradually affects the whole structure of the vortex 
lattice. It affects the form factor of /i(r) and the magnetic 
field distribution function. The former is detected by the 
SANS experiment, and the latter corresponds to the res- 
onance line shape in the /xSR and NMR experiments. At 
higher field, we also have to consider the vortex lattice 
effect on the vortex structure. The orientation of the 
vortex lattice (i.e., the relative angle between the a axis 
and the NN vortex direction) is important for the vortex 
structure. 

The vortex core radius decreases similarly for both the 
s-wave and the dx2 -y2-w&Ye pairing with increasing field. 
In the triangular vortex lattice case, the six NN directions 



are not equivalent in the vortex structure for the dx2-y2- 
wave pairing. This may be the origin of the distortion 
from the triangular lattice. We showed the influence of 
the (ia;2_j,2-wave nature and the vortex lattice effect on 
the LDOS, and discussed the field dependence of the spa- 
tially averaged DOS. The spatial structure of the induced 
s-wave and da;y-wave order parameter in the dj.2_j,2-wave 
superconductor was also presented. We note that the ef- 
fect of the induced order parameter was not included in 
our discussion of the dr^2_y2-wme nature in Sees. [II- VI. 
The dx2_y2-wsve nature appears even when the induced 
order parameter of the other symmetry is absent. 

The d2;2_j,2-wave nature of the vortex structure A(r), 
J(r) and /i(r) was also studied in the GL theory or the 
London theory by including nonlocal correction terms. 
We note that the GL Iheory only applies by the ex- 
act derivation near Tcx3 For other cases, we have to 
include nonlocal correction terms. We usually consider 
only the fourth order derivative terms as non-local terms, 
which contribute to the vortex structure in the correc- 
tion of the order ln(rc/T). But, at lower temperature, 
we have to include higher order derivative terms. Their 
contribution is a higher order of ln(Tc/r). Further, as 
we increase the amplitude of the order parameter bel- 
low Hrj2, we also have to consider higher order nonlin- 
ear terms of the order parameter than the usual |ApA 
term. It is difficult to consider all contributions of these 
higher order terms by modifying the GL theory. The 
quasi-classical Eilenberger theory automatically includes 
all contributions of the higher order terms. Neverthe- 
less, the da;2_j,2-wave nature and the vortex lattice effect 
obtained by the quasi-classical calculation can be repro- 
duced qualitatively by the GL theory i^iiih the correction 
of the fourth order derivative terms. El'EZi We note that 
the c?a;2_y2-wave nature in the resonance line shape P{h) 
is also reproducedB^I This means that the extended GL 
theory is still useful as a phenomenological theory for the 
qualitative study of the c?2,2_j^2-wave nature in the vortex 
structure. The validity of the study with the nonlocal 
correction term (i.e., fourth ocdec^erivative term) in the 
GL theory or London theoryo £3 can be confirmed by 
a comparison with the result obtained in a microscopic 
theory such as the quasi-classical theory. But, for a quan- 
titative study, we have to use a microscopic theory. 

In the London theory, the amplitude of the pair po- 
tential is treated as a constant outside the small region 
of the vortex core. However, this is only justifiable if the 
field is low enough. In our calculation where T /Tc = 0.5, 
|A(r)| shows a spatial variation at every region of a unit 
cell for H > 0.1Hc2- We have to carefully use an approx- 
imation of the constant amplitude over a wide range of 
the apphed field. 

We expect that the dx2_y2--wave nature presented in 
this paper will be observed in experiments on the dx2-.y2- 
wave superconductors such as the high _Tc supercon- 
ductors and the organic superconductors.c3 We recom- 
mend that the vortex structure be investigated into the 
higher field range, where the (ia;2_„2-wave nature is em- 
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inent. At lower fields, we also consider the effect of 
the gradual transformation from the triangular to the 
square vortex lattice. Borocarbide superconductors such 
as LuNi2B2C are also good candidates for observing the 
dj.2_j^2-wave nature.B These materials are considered to 
be dx2-y2-wave superconductors or s-wave superconduc- 
tors with a highly anisotropic Fermi surface with four- 
fold symmetry. Even for the latter s-wave case, we ob- 
tain a similar vortex structure to the (1^2 _y2 -wave pair- 
ing case by the k-dependence of vp instead of the k- 
dependence of (/){d)L3 There, the k-dircction of large v-p 
plays the same role as the small |0(^?)| direction. In 
experiments on LuNi2B2C, the fourfold-symmetric be- 
havior of Hc2 is reportecl-when the field direction is ro- 
tated within the ab plane.E^I From the Hc2 anisotropy, the 
magnitude of the dx2 ^y2-wave-like behavior is estimated 
as im'ivk - 3<4))k/4('^W'« + «i,)>|=0.43, 
where (• • means the average on the Fermi surface. 
This magnitude is comparable to 0.25 of the (i^2_j^2-wave 
pairing case (0(6') ~ -\/2cos20 and an isotropic 2D Fermi 
surface). So, the large (i^2_j^2 -wave-like anisotropy effect 
is expected in the vortex structure. The gradual tran- 
sition from the triangular to the sq uare .vortex lattice 
is also observed with increasing field.E§E3 As /xSR and 
SANS experiments have begun on these materials, we 
expect that the c?^2_j^2 -wave-like nature will be detected 
there. 

For the other symmetry of the anisotropic super- 
conductors such as heavy fermion superconductors or 
Sr2Ru04, if the k-dependence of their pair potential or 
the Fermi surface structure has a large anisotropy in the 
directions perpendicular to the applied magnetic field, 
the anisotropy effect appears by the same origin as that 
discussed in this paper. It may give important informa- 
tion for determining the symmetry of the order parame- 
ter. 
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FIG. 2. Spatial variation of the vortex structure at low field H /Hc2~0.021 for the ci45° case. The a axis (i.e., 0° direction) 
is horizontal, (a) Amplitude of the pair potential. The contour lines of |A(r)|/Ao are plotted. From the center, the lines are for 
|A(r)|/Ao=0.1, 0.2,..., 0.9. Contour lines show a fourfold-symmetric shape around each vortex, (b) Amplitude of the current. 
Contour lines of | J(r)|/(c(^o/47r^o) are plotted. There are four small peaks around each vortex, (c) Internal field distribution. 
Contour lines of h{r) / {(j}o / £,q) are plotted, (d)-(f) are, respectively, the profiles of (a)-(c) along the lines OA (solid line), OB 
(dotted line), OC (dot-dashed line) and AC (thick line) of Fig. |[ 
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• • 

FIG. 3. Square vortex lattice configuration in a contour plot for the s-wave and the d45° cases. The a axis (i.e., 0° direction) 
is horizontal. For the dO° case, this configuration is rotated by 45°. The vortex centers are shown by solid circles. A square 
enclosed by dashed lines represents the Wigner-Seitz cell of the vortex lattice. In Figs. ^ ^, ^ and ^, the profile of the vortex 
structure is presented along the lines OA, OB, OC and AC. 
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FIG. 4. Spatial variation of the pair potential at high field H /Hc2~0.54 for the s-wave (a), the d45° (b) and the dO° (c) 
cases. The contour lines of |A(r)|/Ao are plotted. From the center, the lines are for |A(r)|/Ao=0.1, 0.2,..., 0.7. The a axis 
(i.e., 0° direction) is horizontal, (d)-(f) are, respectively, profiles of (a)-(c) along the lines OA, OB, OC and AC of Fig. |^. 




(a) H/Hc, (b) H/Hc2 

FIG. 5. Field dependence of Ann and Annn, which are normalized by Aq. We define Ann (Annn) as the maximum of 
|A(r)| along the NN (NNN) direction, (a) s-wave case, (b) ci45° and dO° cases. 
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FIG. 6. Field dependence of the vortex core radius (a) and ^2 (b) for the s-wave (o), the d45° (o) and the dQ° 



(□) cases 

in the square lattice. They are normalized by ^o- The radius ^1 is defined from the initial slope of the pair potential by setting 
as |A(r)| = AnnT"/^!. The radius ^2 is defined from the maximum of the screening current |J(r)|. Both and $2 show a 
similar field dependence. The triangular lattice case of the s-wave pairing is also plotted (•). Its field dependence is almost the 
same as the square lattice case of the s-wave pairing. 
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FIG. 7. Spatial variation of the internal current at high field H/Hc2=0.54: for the s-wavc (a), the d45° (b) and the dO° (c) 
cases. Contour lines of |J(r)|/(c0o/47r5o) are plotted, (d)-(f) are, respectively, profiles of (a)-(c). 
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FIG. 9. Field dependence of the dominant form factor hifi for the s-wave (o), the d45° (o) and the dO° (□) cases in the 
square lattice. The triangular lattice case of the s-wave pairing is also plotted (•). 
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FIG. 10. Field dependence ol higher order form factors. We plot /i2,o//ii,o (a), /ii.i/fti.o (b), fta.o/Zii.o (c) and /i2,i//u,() (d) 
for the s-wave (o), the d45° (o) and the dO° (□) cases in the square lattice. The triangular lattice case of the s-wave pairing is 
also plotted for reference (•). 
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FIG. 11. Magnetic field distribution function P{h) of the square vortex lattice case at H/Hc2=0.021 (a), 0.15 (b) and 0.54 
(c) for the s-wave, the d45° and the dO° cases, h is normalized by 4>o/&- In (a); the lines for the triangular vortex lattice are 
also plotted. 



15 





0.5 1 0.5 1 

(a) H/Hc, (b) H/Hc2 

FIG. 13. Field dependence of variance cr (a) and skewness parameter a (b) for the s-wave (o), the d45° (<>) and the dO° (□) 
cases in the square lattice. The triangular lattice case of the s-wave pairing is also plotted (•). 
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(cl) (ft) (f) 

FIG. 14. Contour plot of the vortex structure in the triangular vortex lattice for the d^2_j,2-wave pairing. As for the 
orientation, Oq = 0° [(a)-(c)] and 45° [(d)-(f)]. The amplitude |A(r)| [(a) and (d)], the current |J(r)l [(b) and (e)], and the 
internal field h{r) [(c) and (f)] at H/Hc2 ~ 0.15 are presented. The a axis (i.e., 0° direction) is horizontal. To show the 
distortion of the contour lines clearly in (a) and (d), we also present the contour line for |A(r)|/Ao = 0.89 in addition to the 
lines of 0.1, • • ■, 0.9. 
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FIG. 15. Magnetic field distribution function P{h) of the triangular vortex lattice case at H/Hc2~ 0.15 for the s-wave, d45° 
and dO° cases. 




(ii) ih) (C) 

FIG. 16. Spatial variation of the LDOS, N{E,r)/No, for the s-wave pairing at H/H^2= 0.15. E/Ao = (a), 0.4 (b), and 0.8 
(c). The region of the Wigner-Seitz cell shown in Fig. His presented. In these figures, the peak is truncated at N{E, r)/No — 2. 
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FIG. 18. Spectrum of the spatially averaged DOS, N{E)/No, at /////c2=0.021, 0.15, 0.54 and 0.75 for the s-wave (a) and 
the d^2 _y2-wa,ve (b) pairings. 




FIG. 19. Spatial variation of the induced s-wave and d^^y-wave components at low field H/Hc2 = 0.021 [(a) and (b)] and 
high field H/H,2 = 0.54 [(c) and (d)]. The amplitude |A,(r)l/lV;/'l/ 2.^,2 1 [(a) and (c)] and \A^y{r)\/\V^y/V^2_y2\ [(b) and (d)] 
are presented in the Wigner-Seitz cell region shown in Fig. pi 
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FIG. 20. Spatial structure of the relative phase for the induced s-wave component [(a) and (c)] and the dajy-wave component 
[(b) and (d)] at low field [(a) and (b)] and high field [(c) and (d)]. The position of the singularity and its winding number are 
presented schematically for arg(As(r)/A^2_j^2(r)) and arg(A2;j,(r)/A^2_j,2(r)). • and o mean the singularities with winding 
numbers +1 and —1, respectively. The solid lines show the cut of the phase, where the phase jumps from — tt to n. The dotted 
line shows the Wigner-Seitz cell shown in Fig. 0. 
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